Exact self-similar solutions in Born— Infeld theory 
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We present a new class of exact self-similar solutions possessing cylindrical or spherical symmetry 
in Born-Infeld theory. A cylindrically symmetric solution describes the propagation of a cylindrical 
electromagnetic disturbance in a constant background magnetic field in Born-Infeld electrodynam- 
ics. We show that this solution corresponds to vacuum breakdown and the subsequent propagation 
of an electron-positron avalanche. The proposed method of finding exact analytical solutions can be 
generalized to the model of a spherically symmetric scalar Born-Infeld field in the (n+l)-dimensional 
Minkowski space-time. As an example, the case n = 3 is discussed. 
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Born-Infeld electrodynamics was proposed in the 
1930s as a nonlinear generalization of Maxwell electro- 
magnetism [l[. The central idea of Born and Infeld was 
to create a relativistic field theory which could admit a 
finite energy classical solution describing an elementary 
electric charge. Since the appearance of quantum electro- 
dynamics (QED), the interest in classical field theory has 
faded considerably. However, in that time, Born-Infeld 
(BI) theory received a fairly unexpected development. It 
was related to Heisenberg's suggestion to describe multi- 
ple 7r-meson emission in high-energy hadronic collisions 
as an expansion of a nonlinear wave packet in the con- 
text of the effective scalar BI field [2]. This classical 
Heisenberg's model has proven to be very fruitful and 
continues to be developed in the physics of quark-gluon 
plasma US]. 

In the past decades, since its rediscovery in low energy 
limit of string theories [B-floj]. BI electrodynamics has 
been studied extensively by a large number of workers 
(see, e.g., [sj-TToj and references therein). A substantial 
degree of interest in the subject has also been stimulated 
by recent advances in laser technology. Modern high- 
power lasers make it possible to reach an intensity level 
of 10 26 -10 28 W/cm 2 . The corresponding field strength 
is sufficient to make nonlinear electrodynamic effects in 
vacuum measurable The experimental verifica- 



tion of BI electrodynamics could serve as an important 
evidence in favor of string theory, and various ways for 
laboratory testing of the corresponding nonlinear effects 
are currently discussed in the literature 19l-23|. 

It is known that BI electrodynamics has unique prop- 
erties (causal propagation and the absence of birefrin- 
gence 24J-26}) among other relativistic nonlinear theo- 
ries of the electromagnetic field. Born-Infeld theory pos- 
sesses an impressive mathematical beauty and has al- 
ready found numerous applications in different branches 
of physics. 

Because of a rather complicated nonlinearity of the BI 
field equations, only several exact solutions are known in 
this theory. These are the point charge solution found 
by the creators of the theory [lj], 2D electrostatic solu- 



tions [27|, |2j| , and plane wave solutions [29h32| . The 
problem of finding new exact solutions in BI electrody- 
namics is very topical. It is well known that BI equations 
admit the existence of exact static singular solutions with 
finite total energy (the so-called Blon solutions). The 
simplest cases are cylindrically or spherically symmetric 
solutions with a singularity on the axis or at the origin, 
respectively. Recently, such solutions has received much 
attention in string or M theory @, H| . It should be noted, 
however, that the properties of nonstationary solutions, 
which describe the propagation of cylindrical and spheri- 
cal wave disturbances, remain poorly studied not only in 
BI theory, but also in Maxwell electrodynamics of non- 
linear media. Some exact axisymmetric solutions of the 
Maxwell equations in a nonlinear medium have recently 
been found in 33-35j|. In this work, we obtain exact 



self-similar solutions possessing cylindrical or spherical 
symmetry in BI theory. 

The Lagrangian density of BI electrodynamics is given 

by 



(Air)-^ 2 (l - v 7 ! - b 



- 2 I - b~ 4 J 2 



(1) 



where b is Born's constant I 



\F^F^ LU = E 2 -B 2 
and J = — \F iiV J :ilv = E • B are the Poincare invariants, 
and F^ u and are the electromagnetic-field tensor 

and the dual tensor, respectively. The electromagnetic 
field equations following from Lagrangian density (Q~]) for- 
mally coincide with the Maxwell equations if the con- 
stitutive relations have the form D = 4irdL/d~E and 
H = — 4ir dL/dti. We introduce a cylindrical coordinate 
system (r, <fr, z) and assume that the fields are indepen- 
dent of <f> a- n d z. Then the BI equations admit solutions 
in which only the E z and B^ components are nonzero. 
Denoting these components as E and B, respectively, and 
using the fact that J vanishes in this case, we can write 
the field equations in the form 



d r H 



V H = d T D, d r E = d T B, 



(2) 



where D = E(l — b^I)- 1 / 2 , H = B(1- ft" 2 /)- 1 / 2 , I = 
E 2 — B 2 , t — ct, and c is the speed of light. The second 



2 



equation in system ([2]) is satisfied by putting 

E=-bd T ip, B = -bd r ip, (3) 

where if) is the normalized (to b) z-component of the vec- 
tor potential. Assuming that I ^ b 2 , from the first equa- 
tion in system ([2]) we obtain 

[1 - (d T i)) 2 ]d^ + 2d r ^d T lPd r (d T ^) - [1 + (d r i)) 2 }d 2 ^ 
+ (n - l)^ 1 [1 - {dri,) 2 + (d r i)) 2 ]d r ^ = 0, (4) 

where n = 2 in the case considered. Here, the integer pa- 
rameter n is introduced for the following reasons. Equa- 
tion (4J can be considered as the Euler-Lagrange equa- 
tion for the scalar field ip{r, t) in the (n + l)-dimensional 
Minkowski space-time, which follows from the action 

S = [[l-(d T ^) 2 + {d r iP) 2 ] 1/2 r n ~ 1 drdT, (5) 



where r — y x \ + • • ■ + x n ■ Although our main atten- 
tion will be focused on the cylindrical symmetry in BI 
electrodynamics (n = 2 and r = \/ x 2 4- y 2 ), we will also 
consider the spherical symmetry which corresponds to 
the case where n — 3 and r = \/ x 2 + y 2 + z 1 . The gen- 
eralization to the higher dimensions is straightforward. 
Equation (4} admits self-similar solutions of the form 



tp = ru(s), s = rr 1 . 



(6) 



Substituting Eq. (J6j) into Eq. ((4]) yields the ordinary dif- 
ferential equation 

{s 2 -u 2 -l)u" + (n-l)(u-su')[l-(u') 2 + {u-su') 2 ] = 0, 

(7) 

where the prime denotes the derivative with respect to s. 
We will seek an exact solution of Eq. (7]) in parametric 
form: 



£ 1 / 2 cosh a, s = £ 1//2 sinha. 



(8) 



Here, a — J <l>(£)e?£ + q, where q is an arbitrary integra- 
tion constant. Substituting expressions (SJ into Eq. ([7]) 
and using the formulas 

< = 4/4, < s = (44' £ -44«)/(4) 3 ; (9) 

we arrive at the Bernoulli equation 

d£, ^ LV ^ J 2f(f+l) v ' 



(11) 



Integration of Eq. (10]) gives 

1 (C + l) 7 



$ = ± 



where 7 = (n - l)/2, Xn (f ) = (4 + 1)" - (p + n)£ is a 
polynomial of order n, and p is an integration constant. 



Restricting ourselves to consideration only of the simplest 
cases n = 2 and n = 3, we write down 



and 



X 3 (£)=f + 3£ 2 -K + l. 



(12) 



(13) 



Thus, Eqs. (jH), ((8|), and (fTT|) give an exact solution of 
Eq. Q. From these expressions, we have 



d T 1p = ± 

= ± 



\/Xn(£) cosh a + (£ + l) 7 sinha 
\J Xn (0 sinh a + (£ + l) 7 cosh a' 

VXn(0 sinha + (£ + l) 7 cosh a 



(14) 



Now we should examine what physically meaningful so- 
lutions can be obtained by appropriately choosing the 
arbitrary constants p and q. 

Cylindrical symmetry. Let us consider the following 
representation of the quantity a in the case n = 2: 



1 /■« VF+Tdf 

a = ±- 



(15) 



where 6 = p/2 - y/p 2 /4-l, 6 - p/2 + y/p 2 /4-l, 
p > 2, and £ > £2- Note that £1 and £2 are the real- valued 
roots of the polynomial x%(0 m ^q. (|12p . Reduction of 
the elliptic integral in Eq. (IB"]) to the standard form gives 

« = ± J-^r Ui - 6)n(/3, i/, k) + (6 + i)F08, k)] , 

(16) 

where F and n are incomplete elliptic integrals of the first 
and third kinds, respectively, f3 — — £2)/(£ — £1), 
v = £1/62 , and fc = V(6 + l)/(6 + !)• We use the 
following notation for the functions n and F: 



n(/9,i/,fc) 



/3 



(1 - k 2 ) V(i - c 2 )(i - w 



(17) 



and F(f3,k) = H(/3,0,k). The components of the ax- 
isymmetric electromagnetic field in BI electrodynamics 
are given by Eqs. (3), (12]), (14]), and (16]). With the 
identity ct/r = i; 1 / 2 sinha, these formulas determine J5 
and _B as functions of the radial coordinate and time via 
the parameter £ (£2 < £ < 00). 

The energy conservation law 9 t W + V • S = 0, with 
the energy density 



b 2 fl + b- 2 B 2 
W = — , - 1 



4vr VVl - 6 -2 / 
and the Poynting vector 

„ c EB 



4tt VI -b- 2 r 



(18) 



(19) 
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FIG. 1. Normalized field components _E = _E/b, B = B/6, 
D = -D/&, and = _H~/6 as functions of r/r for p = 2.2 (a,c) 
and p = 6 (b,d). 
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FIG. 2. Normalized energy density, the radial component of 
the Poynting vector, and the invariant I (AirW/b 2 , 47r£/c6 2 , 
and I/b 2 , respectively) as functions of r/r for p — 2.2 (a) and 
p = 6 (b). It is seen that I/b 2 —¥ 1 for r/r — > ±oo. 



can easily be derived directly from the field equations © . 

Figure 1 shows the results of calculations of E and B 
by Eqs. © and CEU) for p = 2.2 and p = 6. The quan- 
tities -D and are also presented in this figure for the 
same values of p. The plots of Fig. 1 can be considered 
as oscillograms of the field quantities at a fixed point 
r = const =/= 0. The branch of f 1 / 2 in Eq. flSJ and the 
signs in Eq. ()14|) were chosen to provide the continuity 
of E and B at r = for any r / and ensure that 
the radial component E of the Poynting vector is nega- 
tive for r < and positive for r > 0. Figure 2 shows 
the energy density W, the radial component S of the 
Poynting vector, and the invariant / as functions of r/r. 
Figure 3 presents the snapshots of E and B as functions 
of r for p = 2.2 and various values of r. Since there is no 
characteristic spatial scale in the problem considered, the 
radial coordinate is given in arbitrary units. It is seen in 




10 20 30 

r (arb. units) 



r (arb. units) 



FIG. 3. Normalized electric and magnetic fields (E = E/b 
and B — B/b, respectively) as functions of r for p = 2.2 and 
various values of r: n < and ri < T2 < (a) and T3 > 
and T4 > tz (b). 



Figs. 1 and 3 that the field quantities described by the 
obtained solution are single- valued continuous functions 
at any space-time point, except for the points on the sym- 
metry axis. The solution describes the propagation of a 
cylindrically symmetric disturbance in the azimuthally 
magnetized vacuum. For r < 0, the cylindrical electro- 
magnetic wave converges to the axis and the wave profile 
becomes steeper [see Fig. 3(a)]. As a result, a shock wave 
forms at r = and r = 0. At this time instant, the elec- 
tric field and the energy-flow direction reverse their signs, 
so that for r > we observe propagation of a divergent 
cylindrical wave [see Fig. 3(b)]. Although the fields E 
and B are everywhere finite, it can be shown that D, 
H, and W behave as const x r _1 for r — > and t^0. 
Such a singularity is responsible for a linear increase in 
the dependences D(r/r), H(r/r), W(r/r), and E(r/r) 
in Figs. 1(c), 1(d), 2(a), and 2(b). It is also seen in these 
figures that (i) an increase in W corresponds to the pos- 
itive values of the invariant /, (ii) the quantities D, H, 
and W diverge not only at r = const =^ and r — > 0, 
but also at r = const =^ and r — > 00, and (iii) the lo- 
cal group velocity v g = |S|/W does not exceed c. The 
energy density is everywhere integrable. 

The presence of the above singularity allows us to pro- 
pose a physical interpretation of the obtained exact so- 
lution as that due to a distributional source on the axis, 
which is similar to static Blon solutions (I,®. The di- 
vergent cylindrical wave [see Fig. 3(b)] can be excited by 
"switching-on" of a delta-function source on the axis at 
the time instant r = 0. Since charges and currents can 
not be specified independently of the field in nonlinear 
BI electrodynamics [l(, the existence of such a source is 
inseparably related to the presence of a constant back- 
ground azimuthal magnetic field B{r) = B, m = b^f& > b 
at r = and r / 0. On the axis, the chosen source 
supports constant fields E(r — 0, r > 0) = — Eq and 
B(r = 0, r > 0) = B such that I = E% - B$ = b 2 
and (-D, H)—> 00. This limiting field state existing in BI 
electrodynamics has received an interpretation in string 
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theory as a divergence in the rate of pair production of 
open strings 0, [36| . An analogous effect follows from the 
QED model of vacuum polarization. In QED, the elec- 
tric field the strength of which is close to the Schwinger 
limit may cause electron-positron pair creation from vac- 
uum [37J . The necessary condition for the Schwinger pair 
creation process is / > 17], [37| . The interaction of the 



created electrons and positrons with a sufficiently strong 
field can lead to production of multiple new particles and 
avalanche-like vacuum breakdown (Hi, HI]. Thus, we can 
state that the obtained solution describes radial expan- 
sion of the electron-positron plasma bunch. Due to an 
avalanche-like electromagnetic cascade, we have infinite 
polarization and magnetization of vacuum. This process, 
which is created on the axis, progressively propagates 
in the whole space. From the classical viewpoint, such 
an expansion can be explained intuitively as a drift of 
charged particles in the crossed fields E z and B$. 

Spherical symmetry. Since continuous tangential vec- 
tor field on a sphere cannot depend only on the radial 
coordinate, this case obviously has no direct bearing on 
electrodynamics. However, the model of a spherically 
symmetric scalar field with BI action ([5]) is studied in- 
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tensely in connection with string/M theory [4C 

In what follows, we consider the solution of Eq. Q 
with n = 3 in a way similar to that used for the axial 
symmetry. For n = 3, the constants p and q can be 
chosen so that the quantity a is given by 



f 



±- 



£n/(£-6)(£-6)(£-6)' 



(20) 



where £i, £2, and £3 are the real- valued roots of the 
polynomial X3(0 m Eq. ([TB")) such that £3 > £2 > £i- 
Since complete analysis of the cubic equation is cumber- 
some, we consider only a single value p = 5 in Eq. (|13p. 
In this case, the above-mentioned roots are equal to 
£1 = — 2 — V5) £2 = — 2 + V5, and £3 = 1, so that a 
can be represented as 



±- 



[(&-i)n(/3,6,A0 + fo 



l)F(/3,k)}, 

(21) 

where /3 = !)/(£- 6) and k = ^(6 -&)/(! - 

The quantities d T ip and d r ip are given by Eqs. (Ti"3"l) . (Tl"4"l) . 
and (|2T1) . The energy density and flux can readily be 
found from Eqs. (|18p and (Tl9|) by representing E and 
B in terms of the derivatives of ip in accordance with 
Eq. (j3|). The results of calculations of d T ip, d r ip, W, and 
E are shown in Fig. 4. It is seen in the figure that by 
analogy with the cylindrical case, the solution describes 
the propagation of a spherically symmetric disturbance in 
a constant background field. It can be shown that W ~ 
const x r~ 2 for r — > (r ^ 0) and, hence, the singularity 
at the origin is integrable. The physical constraint v g = 
\E\/W < c is also satisfied. 
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FIG. 4. Quantities <9 t i/j and (9 r f/> of the spherically symmetric 
BI field (a) and the energy density and the radial component 
of the energy flux (b) as functions of r/r. 



In conclusion, we emphasize that the obtained exact 
solutions exist essentially due to BI nonlinearity and spa- 
tial symmetry. Maxwell equations (b — > 00) as well as 
BI equations in flat geometry [n = 1 in Eq. @] obvi- 
ously do not allow the existence of such solutions. Al- 
though the obtained solutions do not have a finite total 
energy, their energy density is locally integrable. Be- 
cause of this, by suitably cutting off, the obtained so- 
lutions, as, e.g., plane wave solutions, may provide use- 
ful approximations to solutions with finite total energy. 
Moreover, BI wave propagation in a constant background 
field represents considerable interest in string theory [Tol j . 
Finally, we note that Eq. (Q} also admits self-similar solu- 
tions i[> = tu(t/t). The parametrization u = tanha and 
rr _1 = (-£) _1 / 2 (cosha) _1 with a = J$(£)d£ + q leads 
again to Eq. (|I0[) . Because of discontinuity or ambiguity, 
we failed to give some of these solutions any physical in- 
terpretation. However, it can be assumed that among the 
whole set of partial solutions, there may exist physically 
meaningful ones. Since they are determined by three gov- 
erning parameters (n and two integration constants) of 
the problem, an appropriate choice of these parameters 
needs further studies. 
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